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Phenomenological approach is developed in the theory of spin-valve type ferromagnetic junctions 
to describe exchange switching by current flowing perpendicular to interfaces. Forward and backward 
current switching effects are described and they may be principally different in nature. Mobile 
electron spins are considered as free in all the contacting ferromagnetic layers. Joint action of the 
following two current effects is investigated: the nonequilibrium longitudinal spin-injection effective 
field and the transverse spin-transfer surface torque. New vector boundary conditions are derived, 
which represent, in essence, the continuity condition for different spin fluxes. Dispersion relation for 
fluctuations is derived and solved for a junction model having spatially localized spin transfer torque: 
depth of the torque penetration into the free layer is assumed much smaller than the total free layer 
thickness. Some critical value kq of the well known Gilbert damping constant k is established for the 
first time. Spin transfer torque dominates in the instability threshold determination for small enough 
constants k < kq, while the spin-injection effective field dominates for k > kq. Typical estimation 
gives Ko ~ 3 X 10~^ and it shows that the both mechanisms mentioned may be responsible for the 
observable threshold in ferromagnetic metal free layer used at room temperature. Fine interplay 
between spin transfer torque and spin injection is necessary to provide a hysteretic behavior of the 
resistance versus current dependence. The state diagrams building up show the possibility of non- 
stationary (time dependent) nonlinear states arising due to instability development. Calculations 
lead to the instability increment values of the order of < 10^° s~^ and they are really containing 
the contributions from the two mechanisms mentioned. Spin wave resonance frequency spectrum 
softening occurs under the current growing to the instability threshold. Magnetization fiuctuations 
above the threshold rise oscillatory with time for small damping (k < kq), but rise aperiodically 
and much more rapid for large damping {k > kq). 

PACS numbers: 72.25.Ba, 72.25.Hg, 75.47.-m 



I. INTRODUCTION 

Spin valve type magnetic junctions are of very interest 
for investigations and subsequent applications. A num- 
ber of significant effects were predicted and experimen- 
tally revealed in the junctions last years. Magnetization 
instability was predicted in Refs. The instability 

is driven by the current perpendicular to plane (CPP 
regime) due to near interface transformation of a trans- 
verse part of mobile electron spin flux to the lattice mag- 
netization flux. This mechanism is usually referred to 
as a surface transverse spin transfer (TST) torque. An- 
other mechanism was proposed and theory elaborated in 
Refs. 3 4,5 6. This second mechanism is due to nonequi- 
librium longitudinal spin injection (LSI) by current into 
the bulk of a ferromagnetic layer. The injected spins 
create current dependent sd exchange in nature effective 
field, which acts on the lattice magnetization and may 
stimulate reorientation phase transition. Theory of the 
joint action of the two mechanisms mentioned was devel- 
oped recently^ (see also preprintaSiS, and Ref . ITo|) 

Current driven instability was revealed experimen- 
tally, being manifested in resistance jumps'^^ and in mi- 
crowave emission. Many interesting subsequent exper- 
iments confirmed the first observations in principle and 
gave a lot of additional significant information (see, e.g., 
Refs. Il3lll4|[l5il6»17il8^) . The most of experiments con- 



sidered trilayers Co/Cu/Co as a constitutive working 
part of the junction. The first Co layer has usually pinned 
lattice due to large its own thickness or existence of addi- 
tional pinning layers. The second Co layer is taken thin 
enough to have free lattice magnetization, which may be 
controlled by external field or current. It is significant to 
remark that thickness L of the second free layer cannot 
be done less than ~ 2 nm to prevent a non-continuous 
structure appearing. Typically, the thicknesses in ex- 
periments were chosen within the range L ^ 2.5-10 nm. 

Success in experimental study and in applications 
stimulated an activity to achieve more detailed under- 
standing of the fundamental idea of TST process intro- 
duced in Refs. Ill2l Contrary to this, in the present pa- 
per we accept this idea without any additional proof: 
an abrupt transformation appears of transverse mobile 
electron spin flux to the lattice magnetization flux at 
the interface between the layers. Great attention was 
paid already to justify this idea more rigorously start- 
ing f rom the microsco py theory approach (see, e.g., 
Refs. Il9l20l2ll22l23l24j) . However, the phenomenologi- 
cal approach is necessary also because of its more wide 
sight. 

Exactly such an approach we laid in the basis of the 
present paper. In other words, we accept the mobile elec- 
trons cannot perform long time and long distance preces- 
sion around the lattice magnetization because of phase 
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mismatch of individual electrons due to their statistical 
velocity spreadJ^^ As a result, the transverse compo- 
nent of the mobile electron magnetization vector decays 
rapidly and transfers to the lattice of ferromagnet. This 
picture is taken as the base of our consideration and then 
all the phenomenology is built up. 

Such an approach allows us to integrate in the the- 
ory the second mechanism of current action that is LSI. 
Therefore the LSI effective field and TST torque are 
treated simultaneously as the two manifestations of the 
only sd exchange interaction. Additionally, we shall show 
such significant effects as the backward current switching 
and the hysteretic resistance versus current dependence 
may be understood as a fine interplay of the two mecha- 
nisms mentioned. As far as we know, similar unified de- 
scription including the both mechanisms, TST and LSI, 
could not developed yet in the frame of a microscopic 
theory. Moreover, we calculate in the present paper not 
only the threshold current (that was done in many pa- 
pers) but the property of the various magnetic config- 
urations: what configurations remain stable and what 
ones become unstable depending on the current. Sim- 
ilar experimental data were presented in the literature 
before and we have now a possibility to compare some of 
them with the theory. We show (similar to Refs. 7 8) the 
spin wave spectrum of the junction layer should soften 
to zero under the current is growing and tends to the in- 
stability threshold. That is decisive manifestation of the 
LSI mechanism, which specifies it in comparison with the 
TST one. We estimate also the rise time for the unsta- 
ble fluctuations. The minimal rise time may be less than 
~ 0.1ns but only if the LSI mechanism of instability is 
dominated. 

Let us make some more preliminary remarks. The 
first remark concerns to the model of the system under 
consideration. According to the original estimations^'^ 
and subsequent theory calculationSfi2i2Si2L2^ the depth 
of transverse spin flux penetration into the free layer is 
of the order of electron quantum wavelength at the Fermi 
surface Ai? < 1-3 nm and is strongly nonuniform within 
this interval. The spin transfer torque acts, as it is well 
known^i^ in that interval only. Therefore the uniform 
approximation of the torque throughout the whole thick- 
ness L seems to be very approximate at least for large 
enough L > 4-10 nm. Because of the model of the torque, 
which is uniformly spread in the free layer, is traditional 
one, it needs, at least, to be justified. RecentljiSiS we 
proposed another model for the calculations, namely, to 
consider spin transfer torque as a boundary condition at 
the pinned/free layer interface. In essence, this model 
is an opposite limiting case with respect to the uniform 
torque model. The obvious advantage of our model is 
in its complete non-sensitivity to details of the processes 
inside the sublayer adjusted to the interface and having 
thickness Xp. More detailed comparison of these two 
models will be discussed in Appendix. 

The second remark concerns to the external magnetic 
field infiuence. We assume the junction may be placed in 



some external magnetic field H, which lies in the plane 
of junction interfaces. But this field is supposed to be 
less in magnitude in comparison with the layer magnetic 
anisotropy field Ha, that is |H| = < Ha- Such a sup- 
position is accepted for simplicity because we want focus 
our attention here to current (not to field) driven mag- 
netization reversal processes. Recent paper32Si2& show at 
large enough fields H > AnM ^ Ha the situation may 
become much more complicated. 



II. MODEL CONSIDERED 

We will use further the most elementary model of a 
ferromagnetic metal based on the well known concept of 
free s-electron (quasiparticle) ideal gas interacting with 
the lattice"^^ bound d-electron magnetization due to the 
so called sd exchange interactioni^ Some more advanced 
(and more complicated!) models may be used, of course. 
However, on our opinion, the most actual problem now 
is to explain a number of experimental facts concerning 
to switching of ferromagnetic nano junctions. We believe 
it may be done in the frame of the simple model^S 

We assume the junction consists of two ferromagnetic 
metal layers and one very thin nonmagnetic spacer in 
between (shown as a heavy line in Fig. ^ . To close the 
electric circuit, a nonmagnetic metal layer 3 is necessary. 
The ferromagnetic metal layer 1 has pinned orientation 
of the lattice magnetization Mi, but the mobile electron 
magnetization mi is, in general, non-pinned. Another 
ferromagnetic metal layer 2 has free lattice magnetization 
M2 and mobile electron magnetization m2, so that the 
magnetizations direction can be changed by an external 
magnetic field H or spin-polarized current density j. 

Plane a; = is the interface between the layers 1 and 2 
and lattice magnetization Mi lies within the plane. Vec- 
tors Ml and M2 make an angle x of arbitrary value 
under the current is absent, that is under j = 0. Af- 
ter the current is turned on, the electrons flow from the 
layer 1 into the layer 2 or vise versa and then appear 
in a non-stationary quantum state and "walk" between 
the spin energy subbands of the corresponding layer. For 
forward current, when the inequality j /e > is valid, 
vectors m.2 and M2 perform a precession (see Fig.QJ. As 
it was shown first in Refs. QQ, angle of the precession 
decreases with coordinate x increasing and tends to zero 
at a; « Af, where Xp is an electron quantum wavelength 
at the Fermi surface. This is because of electron velocity 
statistical spread and no relaxation processes are needed 
to provide such a behavior. More detailed discussion may 
be found also in the recent preprinti^ Analogous "preces- 
sion" region appears, of course, near the interface a: = 
for backward electron current, that is for j/e < 0. How- 
ever, the region will be localized not in the layer 1 but 
again in the region < x < Af of the layer 2. It appears 
due to the complete pinning of the lattice in the layer 
1. This feature is better to discuss a little below in the 
paper. 
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FIG. 1: (Color online) A scheme of magnetic junction to illustrate processes in layer 2. The contacting layers are labelled by 
1, 2, 3. The arrows show directions of the following vectors: Mi and mi are magnetizations in the layer 1, M2 and m2 are 
magnetizations in the layer 2, external magnetic field H lies in the junction plane x = 0, and current density is j/e (e is electron 
charge). The vertical dashed lines separate two ranges in the layer 2. In < x < Af range, precession takes place, which is 
shown with ovals. Vector m2 has both longitudinal and transverse components m2\\ and m2±, respectively. The precession 
angle decreases with x increasing. In x > Af range, the precession vanishes, so that only the longitudinal component remains. 



The region < x < Xp was introduced firstly by Slon- 
czewskii and Berger^ and further we will refer to it as 
"SB layer" . Their considerations of the layer were based 
on the assumption that a ballistic transport occurs inside. 
An opposite limiting case, when diffusion transport dom- 
inates, was considered in Refs. 19,23. We follow, in this 
point, the approach of the original worksii*^ The validity 
criterion for this assumption may be written as Xp < Ip, 
where momentum mean free pass length Ip may be esti- 
mated as ~ 1-10 nm for metals at room temperature. As 
for the region x > Xp in the layer 2, the length Ip appears 
indeed very small one. Based on this reason, we will con- 
sider further the electron transport outside the SB layer 
(in the depth of the layer 2) using the diffusion-drift 
limit. 



is effective field, Hq = /9(Mn)n is anisotropy field, (3 be- 
ing dimensionless anisotropy constant (typically, /3 ~ 0.1 
for our magnets) and n being a unit vector along the 
anisotropy axis direction, A is the intralattice inhomo- 
geneous exchange constant (magnetic "stiffness" of the 
lattice), is demagnetization field, Hsd is sd exchange 
effective field. The latter takes the form2& 



VLsd (x) 



SM (x) ' 



(3) 



where 6 {. . .) /(5M (x, t) is a variational derivative, and 
Usd is sd exchange energy, 



III. LATTICE EQUATIONS OF MOTION 

Quantities that have no number indexes (1, 2 or 3) 
we refer further to any of the junction layers. Motion 
of vector M will be described by means of the Landau- 
Lifshitz-Gilbert (LLG) continuity equation: 



— +7[M,H.,,]-- 



M, 



dM 
'dt 



0, 



(1) 



where t is time, 7 is the gyromagnetic ratio, k is dimen- 
sionless lattice damping constant (0 < k ^ 1), 



H 



e// 



H + H„ 



A 



dx"^ 



■ Hrf -I- Hs 



(2) 



U 

Usd = - ai y" mi {x') Ml dx' 
—00 

L 

- a2 y m2 (x') M2 (x') dx', (4) 



parameters ai and 012 being dimensionless sd exchange 
constants (typical estimations are ai ~ 0:2 ~ 10'*-10^). 
Due to the last term in Eq. (j^J, the motion of vectors M2, 
mi and m2 appears to be coupled. Remember, vector 
Ml is pinned completely in our model and cannot be 
moved. 
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IV. SPATIAL DISTRIBUTION OF MOBILE 
ELECTRON MAGNETIZATION 

We describe the motion of vector m by means of con- 
tinuity equation 



dm 
'dt 



d3_ 

dx 



70; [m, M] 



0, 



(5) 



where r is a time of relaxation to the local equihbrium 
value m = toM, M = M/A/ is the unit vector, J is a 
mobile electron magnetization flux density. We do not 
consider here the almost ballistic quantum mechanical 
motion of the mobile electrons inside the SB layer. In- 
stead, we consider the influence of this small thickness 
layer on the motion outside it by means of an appropri- 
ate boundary condition (see below). Diffusion-drift ap- 
proximation is applicable in the regions outside SB layer, 
that is for x > Xp. The equation |(SJ) may be simplified 
significantly in this region. 

First, only longitudinal part m|| = toM may be taken 
into account. Therefore an effective frequency lo of the 
motion is determined by the precession in relatively small 
fields H, Hd, (3M < -/aM = ujsd 3 x 10^^ We 
assume the condition 



(6) 



is valid, which allows us to neglect time derivative in Eq.jS] 
in comparison with the relaxation term. The remaining 
equation (jSJ may be solved then with respect to Am = 
(m — rh) — AtoM. Let us take the typical parameter 
estimations: ai ~ 0:2 ~ 2 x 10*, M ~ 1 x lO^^ G, r ~ 
3 X 10"^'^ s to justify the following simplifying condition: 



OJsdT - 10^ > 1. 

Then Eq. (O acquires the form 

Am= -tM (M-d3/dx 



(7) 



(8) 



Second, spin flux J can be found explicitly in the region 
X > Xp to substitute it into Eq. ||SJ). Electrons in the 
region occupy energy subbands having their own spin 
magnetic moments parallel to M (|) and antiparallcl to 
M (I). Magnetization of the mobile electrons m contains 
contributions from the both subbands and therefore may 
be presented as 



(9) 



where > is Bohr magneton, n|_| are partial electron 
densities in the spin subbands. The total electron density 
n = + ni does not depend on x and t because of 
local neutrality conditions in metal. The magnetization 
flux density J in Eq. ((SJ also can be related with partial 
electric current densities in each subbands, and j^. 
The relation has the following form: 



Ms 



(iT-Ji)M. 



(10) 



The total current density j = j| -t- does not depend on 
X too, because of one-dimensional geometry of our model. 
We use standard diffusion-drift formula for partial cur- 
rents: 



dn 



dx 



(11) 



where fi^^i are partial mobilities, -D|.x are partial diffu- 
sion coefficients and E is an electric field strength inside 
the layer. Now, we should substitute formula Hll|) in 
Eq. ifTUIl and express J via m, using Eq. Q. The cal- 
culations are direct and were done completely in Ref. |5|. 
The result may be written as follows: 



(12) 



where Q — (cj — (T|) / (cr-f -f- CTj) may be under- 
stood as an equilibrium conductivity spin polariza- 
tion parameter with a^_i = en-f j/if^j, and D = 
{(J-\Di + aiD^) I (ct| + CTj) is the effective spin diffusion 
constant. To obtain Eq. (|12|l . an additional assumption 
should be made,^ namely 



— < 1, 

3b 



(13) 



where j\d = enl/r is a characteristic current density 
in the layer 2. With typical parameter values, n ^ 
lO^^cm-^, Z ~ 3 X 10"^ cm, r - 3 x IQ-^^s, we get 
jo ~ 1-6 X 10^° A/cm^. The condition means the current 
disturbs the subband populations rather slightly, so that 
only a low injection level is realized. As it was mentioned 
above, we are interesting in currents j < 10^-10^ A/cm^, 
which have an order of instability thresholds. Therefore, 
the condition H13|l should be well satisfied in our calcula- 
tions. 

After substituting expression H12|l into Eq. (jSJ we ob- 
tain finally the following equation for m: 



d^m Am 
dx'^ P 



= 0, 



(14) 



where Am = m — m. Since the sd exchange gap is fixed 
in our model, m value does not depend on x and t. The 

spin diffusion length is I = V Dr. 

As it was mentioned at the beginning of this section, 
our above presented derivation may be referred to any 
junction layer. Therefore, Eq. H14I) is valid to describe 
any definite layer if we perform the replacements: m — )■ 

"11,2,3, M Ml, 2, 3, Q (91,2,3, D Dl,2,3, T n,2,3, 

and so on. Moreover, we should take M3 = Mq, 7713 = 0, 
Qs — for nonmagnetic layer 3 (see Ref. |a for more 
details). 

To solve Eq. 1)14(1 . we should derive the appropri- 
ate boundary conditions. We accept two sorts of the 
conditions: 1) the continuity condition of the mobile 
electron longitudinal spin fluxes, penetrating through 
the interfaces and 2) the continuity condition for dif- 
ference between the chemical potentials of electrons in 
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the spin subbands. It means the following equalities 
are to be valid for current flowing in forward direc- 
tion (j/e > 0): Ji(-Af)M2(Af) = J2(AF)M2(Ai.) and 
J2(i - e)M2(L - e) = JaiL + e)M2(L - e), as weU the 
following equalities are to be valid for backward direc- 
tion of current {j/e < 0): J2(A_f)Mi = Ji(-A_f)Mi and 
33{L + e)M2{L -e) = J2(L - e)M2(i - e). It may be 
written in an explicit form using the definition H12|l : 



QiJ - -Di-jr— 

e ox 

Q2J - D2^— 

e ox 



-—\f 



(MiM2(Af)) 



Q2J - D2^— 

e ox 



x—L—e 



^ dx 



x—L+e 



(15) 



for j/e > and 



. n dm2 
— Q2J - D2^— 

e Ox 

Qi3 - Di--— 

e ox 



\ (MiM2(Af)) 

X—\f / ^ ^ 



~ dm3 
-^3 



dx 



x=L+e 



Q2] - D2^— 

e Ox 



x—L — e 



(16) 



for j/e < 0. 

The continuity condition for chemical potential differ- 
ences was considered before us but, apparently, only for 
collinear orientations of magnetizations on the two oppo- 
site sides of the boundary (spacer— or domain wall^°). 
We have written the condition for oblique magnetization 
orientations in Ref. |^ Now, we should return to the 
situation of oblique magnetizations because: 1) the ide- 
alization of pinned mobile electrons is removed and 2) 
vectors Mi and M2 may appear noncoUinear. After cal- 



culations analogous to Ref. [sj, we obtain 

NiAmii-Xp) = N2Am2{XF) (MiM2(Af)) , 
N2Am2iL - e) ^ NsAmsiL + e) (17) 

for j/e > and 

NiAmii-Xp) (MiM2(Af)) = A^2Ato2(Af), 
NsAmsiL + e) = N2Am2iL - e) (18) 

for j/e < 0. Here the following definitions are intro- 
duced: N = (2^^5)7^ [(5t(C))-' + (5i(C))-'], the quan- 
titles gij^i are densities of states for spin up and spin down 
electrons, C is equilibrium chemical potential, which is 
constant across each junction layer. 

The equation ()14|1 may be directly solved with the fol- 
lowing result: 



Ami (a;) = rhi exp [ — 



Am2(a;) — rh2 cosh 
Am3{x) = 7713 exp 



x 
h 
L — X 
h 



l2fTi2 sinh . 



(19) 



where four constants mi, rh2, rh'2, should be found 
from the boundary conditions (|15|I - H18() . Moreover, we 
take Ami 0, Am^ far from the interfaces, that is 
for \x\ — > 00. Then we obtain the distribution of mobile 
electron magnetization across the junction considered. 
After some direct calculations, we obtain the following 
results for the forward current j/e > 0: 



m2 



AiBn2-^{Q2+ Qi (MiM2(Af)) - Q2 (cosh A -h 1/2 sinh A) I 



sinh A -f 1^2 cosh AH (^MiM2(Af)J (cosh X + 1^2 sinh A) 



(20) 



12^2 = MiM2(Af) I m2-fiBn2- 

1^1 ^ ' JD2 



Qi MiM2(Af) -Q2 



(21) 



mi 



N- 



MiM2(Af) m2 



r 



(22) 



Here we introduce dimensionless thickness A = L/I2 Matching parameters are, in essence, the ratios of cor- 

and matching parameters for different interfaces vi = responding spin contact resistances. For small param- 
D2I1N1/ Dil2N2, V2 = 0^12^2/ D2l-iN-i. eters (z^i_2 ^ 1) spin flux penetrates slowly through 
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the corresponding interface, that is the injection weak- 
ens greatly. Similar situation is typical for metal- 
semiconductor boundary (see, e.g., Ref. Is^) . Opposite 
case of large enough parameter {i^i^2 ^ 1) corresponds 
to ready flux penetration through the interface, and the 
moderate penetration corresponds to 1/1^2 ~ 1. 

As it is seen clearly from formulae (|19|I - H22() . the cur- 
rent disturbs spin equilibrium deeply inside the junction 
layers at distances of the order of diffusion lengths ^1,2,3- 
Therefore the nonequilibrium spin injection by current is 
indeed a bulk effect. Moreover, the spin injection^ de- 
pends on the angle between the vectors Mi and M2(Ai?). 
It shows explicitly the mobile electron magnetization 
does depend on the lattice magnetization vector M2 di- 
rection. In principal, such dependence arises because the 
mobile electrons represent much more rapid system than 
the magnetic lattice. Therefore the electrons "see" at 
any time the instantaneous lattice configuration and are 
adapted to it. Formally, it is seen from our transfor- 
mation |(SJ) of the basic equation for electrons. The 
transformation is based on the condition ©. The lat- 
ter condition allows neglecting any inertia in the electron 



following the lattice. 

Note finally according to boundary conditions H15|l - 
()18|l. we may simply make a replacement 

(MiM2(Af)) ^ (MiM2(Af)) ' (23) 

to go from the forward to backward current direction. 
Therefore, formulae ^ and remain valid for 

backward current (j/e < 0) after the replacement H23|l is 
done. 



V. CALCULATION OF CURRENT 
DEPENDENT SD EXCHANGE EFFECTIVE 
FIELD H^d(i, x) 

We should calculate variational derivative ||2J) to find 
the field HsdU, x). First, we calculate the sd ex- 
change energy substituting formulae (|19|I -H22 () into ex- 
pression Q and performing the integration over coordi- 
nate X. Then, we obtain for the forward current (j/e > 0) 



J_ 

302 



sinh A + 1^2 coshA H (^MiM2(Ai;-)j (cosh A + 1^2 sinh A) 



1 -1 



sinh A -I- — (cosh A - 1) fMiM2(AF) 



-— fMiM2(AF) 



Q2 

Qi 



MiM2(Af)) 



- [sinh A -I- i^2(coshA 1)] 
(cosh X + 1^2 sinh A) 



MiM2(Af) 



Q2 



(24) 



where C/f^^'' is an equilibrium part of the energy, which 
depends only on modulus Mi and M2. Parameter b — 
aiTiMi/a2T2M2 describes the contribution of the layer 
1 into the energy (0}. Note, that the exchange energy 
for backward current {j/e < 0) may be easily got from 
expression (I24|l by means of the replacement (|23|l . 

After differentiating (|24|l in accordance with Q, the 
equilibrium part of the effective field is coUinear with M 
and does not contribute to LLG equation Nonequi- 
librium part of the field may be written as 



dUsd ^ (MiM2(A^^)) 



9(MiM2(Af 
= hsd'^ihSix - \f), 



5m2{x) 



(25) 

where explicit expressions for the field hgd depend on re- 
sults of differentiation of the energy being different for 
forward and backward currents. The most general ex- 
pressions may be directly obtained but appear too cum- 
bersome. We show therefore some simple partial cases 
only: 



a) forward current [j/e > 0) and a very thin layer 2 
(A « 1) 



hsd = A'sa2"2Qi- — 

JD2 



v+ (^MiM2(Af))' 



X < \vi 



V- (MiM2(Af))' 



2hv (MiM2(Af)) 



(26) 



where v = v\V2 — D^liNi/ Dily,N^] 
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b) backward current (j'/e < 0) and A <C 1 



and jSJ, sum them and rewrite in the form 



hsd = AiBa2"2Qi- — 

JD2 



1 + MiM2(Af) 



1 - (MiM2(Af) 



26zy MiM2(Ai.) 



(27) 



It is interesting to note that the general expression 125|) 
for effective field is proportional to singular function d{x— 
Xf)- It means the field vanishes everywhere except the 
SB layer near the interface x = taken from the side of 
the layer 2. Therefore the magnitude of the effective field 
is a bulk quantity depending on such bulk parameters as 
L, I2, 712, V- But this field is applied near to interface 
only, namely, to the point x = Xp. 

This feature has a significant consequence. The action 
of the field may be described by two equivalent ways. One 
of them is to consider field l|25() as a singular term in the 
bulk equations of motion LLG for the lattice magneti- 
zation. Exactly such a way was chosen in our previous 
preprint^ and article J But much more simple way is to 
include the field into the boundary condition for the lat- 
tice magnetization at the boundary point x = Xp- This 
last way was proposed in Ref. and we develop it here 
in the next section. 



VI. BOUNDARY CONDITIONS FOR LATTICE 
MAGNETIZATION 



We suggest a model, which is opposite to the model 
of uniformly spread spin transfer torque in the layer 2. 
In our model, the spin transfer torque is localized inside 
the SB layer < a; < A_f and may be described as a 
boundary condition. Such a model corresponds to the 
real situation if only the thickness of SB layer Xp occurs 
the smallest length among the other lengths in the system 
(i.e., L > Xp, \fA > Af, and / > A_f). 

The total magnetization flux in the layers 1 and 3 con- 
sists only of mobile electron fluxes because magnetic lat- 
tice in the layer 1 is pinned, while the lattice in the layer 3 
is nonmagnetic one. Therefore, fluxes in these layers may 
be described by means of expressions of the type (|12|l . 
The situation with the layer 2 occurs significantly more 
complicated: mobile electrons and lattice both contribute 
to the flux. Mobile electrons contribute to longitudinal 
part of the flux, which coincides with Eq. (|12|l . while the 
lattice creates a transverse part of the flux. 

To derive all the spin fluxes, which may flow in our 
junction, it is convenient to return to initial equations (^) 



a(M + m) 
dt 



- 



M, 



dx^ 



-7 [M, H']+7a[m, M] - — 



-7[M, H,d] 



dx 



m — m 



0, 



(28) 



where H' = H + -I- H^;. The summation of the initial 
Eqs. CQl and jSJ is necessary to do in our case because the 
mobile electron flux may transfer to the lattice one and 
these fluxes together determine dynamics of the system. 
The second summand in the left hand side of Eq. H28|l 
may be rewritten as follows: 



7^ 



M. 



&^x 



d_ 

dx 



M, 



dm 

dx 



dJ 



M 



dx 



where 



"Af 



M. 



dM 

dx 



(29) 



(30) 



may be considered as a lattice magnetization flux with 
parameter a — jAM having the sense of a lattice magne- 
tization diffusion constant (see Refs. 7 8 for more details). 

The next (third) summand in Eq. (|28|l may be trans- 
formed using formula (|25|) for H^jj to the following form: 



7 [M, B.sd] 



dJ 



sd 



dx 



where 



Jsdix) = jhsdh M2(Af), Ml e{x - Af 



(31) 



(32) 



represents a very interesting new quantity that may be 
called as an "effective field transverse flux" and 0{x — 
Af) is a well known step function: 6{x — Af) = for 
X < Xp and 9{x — Xp) — 1 for x > Xp. Therefore, 
two lattice magnetization fluxes exist: (|30|l and H32() . the 
last appearing due to sd effective field. All the lattice 
fluxes are perpendicular to the vector M2. The fourth 
summand in Eq. (|28|) originates from mobile electron spin 
flux J, which is given by Eq. H12() . As it is seen, the flux 
J is directed along vector M. 

Let us introduce the integration over the near-interface 
regions — Af < a; < — e and £ < x < Af with infinitesimal 
parameter e +0 (see Fig. |2|for illustration): 

— s Af 

(...)-= J {...)dx and {...)+ = J {...)dx (33) 

— Xp e 

and apply these operations to Eq. (|28|l . We should now 
introduce the conditions, which say explicitly that our 
nonmagnetic spacer layer appears nontransparent for the 
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exact is the model. Formally, conditions (|36|l and 1)3 7|) 
are satisfied well in the limit A^? — > 0. 

Then, we apply 1(23) to using Eqs. 1(231), and 
making all the simplifications. We obtain as a result 



^JM(-Ai.)- J(-AF)+J(e) = 0, 



(38) 



J(Af) - J(e) + 3m{Xf) + JsdiXp) = 0. (39) 

We put Jjif (— Ai?) = in Eq. l(SS|) because of Mi magne- 
tization is pinned and, as we assume, cannot be excited 
by the electric current. The conditions and (|^ are 
valid for all the fluxes at the interface x = 0. Analo- 
gously, we may derive the conditions for the second in- 
terface X ^ L. The result is 



J(i + e) - J(i - e) - 3m{L - e) = 0, 



(40) 



FIG. 2: Range of integration l|88|l : the vertical dotted lines 
sliow limits of integration. 



lattice magnetization excitations and fluxes. These con- 
ditions may be written as 



Jm(-£)=0, Jm(£) = 0, 3sdie) = 0, 



(34) 



the latter condition being satisfied automatically because 
of the definition l|32|l . Moreover, we accept the condition 
that our spacer is completely transparent for mobile elec- 
tron spin flux that is no spin scattering processes occurs 
inside. This last condition may be written as 



J(-£)=J(e). 



(35) 



After the operations (|33|l are applied, a number of sim- 
plifications become possible. In particular, two last sum- 
mands in the left hand side of Eq. H28() may be considered 
as a very small in comparison with the first summand 
(since < k ^ 1, ujt <C 1). These relaxation summands 
should be omitted if we want the flux continuity equa- 
tions may be valid approximately. Further, there exist a 
number of "volume" summands, which are proportional 
to the small length Xf- It is desirable these summands 
to be sufficiently small too, namely, 



A, 



d{M + m) 



dt 

« |J(-Af) 
d{M- 



m 



dt 



- 7 [M, H'] + 7a[m, M] 



- 7 [M, H'] + 7a[m, M] 



< |J(AF)+Jsd(AF)| 



-Xf <x< — e 

(36) 



e<x<\F 

(37) 



The conditions H3t)|) and (|37|l characterize an applicability 
of our model: the better they are satisfied, the more 



where we take into account the fiux J^.^ is constant at a: = 
L and therefore falls out from the continuity condition. 
The conditions (|38|I -I|4U () represent in the most general 
vector form the quest boundary conditions. They allow 
us further to find all the necessary solutions of our initial 
equations of motion |^ and Q . 

In particular, we obtain immediately the longitudinal 
boundary conditions (fT^ and ltTH|) from Eqs. 
To show this for forward current j /e > we should sim- 
ply multiply scalar Eqs. (|2HIl-(Eni) by M2- For j/e < 0, 
we multiply Eq. gOl) by M2 and multiply Eqs. 
and H4U|I by Mi. Note, after the multiplication the prod- 
ucts appear Jj\/(AF)Mi ~ sinx and J sd{XF)M.i ^ sinx, 
which are small at x <C 1 and at |7r — x| ^ 1- We shall 
consider in this paper (see Section IVIII|I only small or 
close to TT angles x- It gives us the reason to neglect 
further the products indicated that are proportional to 
sinx- 

The next step is to extract transverse part of the condi- 
tion 139() . According to Eq. (|38|l . this part for the forward 
current j/e > is induced by the component 



Ji(e)= M2(Af), J(-Af), M2(Af) 



and leads to the relation 



J_L(e) = Jm(Af) + Jsd(AF) 



(41) 



(42) 



For backward current j/e < the other electron flux, 
namely, J(Af), induces transverse (perpendicular to Mi) 
component while the flux J(e) appears purely longitudi- 
nal one being coUinear with Mi. This interesrting pro- 
cess goes due to current dependent sd exchange effective 
field AHsd (see Eq. produces strong and locahzed 

influence on M2 by virtue the LLG equation should be 
satisfied. The field AH^d is directed parallel to Mi and 
therefore any transverse component of the mobile elec- 
tron magnetization Am2 inside the SB layer will decay 
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and transfer to the lattice, the transformation being anal- 
ogous, in principal, to TST idea of the Refs. There- 
fore, we may write the following induced transverse com- 
ponent for backward current j/e < 0: 



and 



Ml, 



J(Af), Ml 



J±(Af) = 

which leads from Eq. H39() to the relation 

J_l(Af) = — Ja/(Af) — Jsd(AF) 



(43) 



(44) 



Transverse part of condition l|4()|l occurs the same for 
both forward and backward current directions, j/e > 
and j /e < 0, namely 



Ja/(L-£) = 0. 



(45) 



Now, we will try to represent conditions (|42l) and H44() 
in the forms, which are more convenient for further cal- 
culations. Let us multiply Eqs. and by M2(Af) 
vectorially and use relations ((201, EU, lEJ 

and H12|l . Then, we get for j/e > 



dM2 



dx 
-P- 



= k- Ml, M2(Af) 

M2(A;^), [Mi, M2(A;.) 



(46) 



where two new parameters k and p appear, which arc 
given by formulae 



— rr — Qi] - Di—^ 

aM2 \ e ox 



Q2J 



~ dAmo 
D2 



and 



dx 



x—Xf. 



for j/e > 0, 
(47) 

for j/e < 
(48) 

(49) 



where h^d should be taken from formula H26() for j/e > 
and from (|27|l for j/e < 0. 

Parameter k originates from the contribution of mo- 
bile electron spin fluxes Jij_(e) (|^ and Jx(Af) 1031, and 
therefore it describes the spin transfer torque action. Pa- 
rameter p originates from Jgdi^p) flux H32I) and describes 
action of the spin injection effective field. Let us substi- 
tute the before found expression Ami (see Eqs. (|19|l - 
()22|l ') for j/e > into H47|) and corresponding expression 
Am2 for j/e < into (|48|l to simplify both the expres- 
sions for k. It gives the following formulae simplified for 
the thin layer 2 that is for A <C 1: 



eaMo 



(MiM2(Af))^ 



for j/e > (50) 



k = 



^"^^2 ^ ■ ^.(MiM2(Af: 



1 



for j/e < 0. (51) 



It is interesting to return now to the case of completely 
pinned layer 1 (that is not only the lattice magnetization 
but mobile electron spins also are pinned). We see from 
the first expression (|47|l the effective diffusion coefficient 
should be taken as I?i 0. Then it follows from the 
parameters I'l and v definitions made after the formu- 
lae (EOl), and 1211): ^ 00, v = V1V2 00. Substi- 
tuting these limiting values to expressions H50f) , H51|l and 
to |(2ni), I2II1 we obtain 



, ^J■BQl j f^Bia2QlT2X 

^ ^ ~T't ' P " 

aA'h e av2 



(52) 



for a completely pinned layer 1 and for 
any sign of the current j/e. Therefore, the sign of 
the parameter k does depend on the sign of the current 
whereas parameter p is positive at any sign of the 
current. We shall see a little later this feature conserves 
in general case of non-pinned mobile electron spins (see 
formulae (|59|l and (jHO))) and it plays significant role in 
determining the instability character. 

It remains now to rewrite the boundary condition H45|) 
for the second interface a; = L in a more explicit 
form. Using the definition (|3()|l and the equality 

(m2 • dM.2/dx) = 0, we find 



9M2 



dx 



= 0. 



= L-6 



(53) 



VII. STATIC STATES 



We start now to find solutions of LLG equation 
satisfying our boundary conditions (|46() and (|53|) . We 
confine ourselves further to the situation when the cur- 
rent density j is considered as a fixed parameter given 
by an external circuit (the so called "current controlled 
regime"). Moreover, we assume an external field H is 
applied along the positive direction of z axis (see Fig. ^ , 
and the anisotropy field Ha is parallel to that axis too 
(remember: H < Ha <C AttM). Then in absence of cur- 
rent, when there is no coupling between the layers in our 
model, the layer 2 magnetization M2 is aligned along z 
axis. Remember the vectors Mi and z make, in gen- 
eral, an angle x the model. Therefore we may write: 
Ml = — y sin X + z cos x, where x, y and z are basis vec- 
tors of our coordinate system. After the current is turned 
on, the direction of the vector M2 may change. Our task 
now is to search the vector M2 change due to the current. 

First, consider the static states possible when the fluc- 
tuations are neglected but the current presents, j/e ^ 
0. We denote the static state magnetization as M2. 
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To find the magnetization we would solve the equation 



M2, He// 



0, which follows from Eq. (P). The last 

equation and the boundary conditions are essentially 
nonlinear and nonuniform. Therefore general problem 
appears rather complicated even for finding the static 
states. 

Fortunately, there exists one simple situation, which 
is at the same time the most actual one, namely, the 
situation with the angle x = or x = tt. In the situation 
all the boundary conditions and the equation of motion 
for M2 are satisfied for 



= M2y - 0, M2z = ±1. 



(54) 



The solution written can be easy justified by substitu- 
tion oP H54|) into the equation and boundary conditions. 
We consider further only these last indicated values of 
the angle x and therefore investigate the stability only 
for the static states ()54(l . For definiteness, we will con- 
sider only the angle x = tt and therefore have 



Ml = -z. 



(55) 



As it is seen from Eqs. and ifS^ . two situations 

are possible: either vectors M2 and Mi are parallel each 
other (P state) or they are antiparallel (AP state). We 
shall consider further the both possibilities. 



VIII. FLUCTUATIONS. DISPERSION 
RELATION 

Passing to consideration of fluctuations AM2, we in- 
troduce them by the following relation: 



M, = ±s 



AM2, 



(56) 



where AM2 = 



AM. 



^ 1. We linearize the LLG 



equation with respect to AM2 and, according to H25|l . 
lay AHsd = outside the SB layer, that is for x > 
Xp. The demagnetization field takes the form = 
— 47rM2AM23;X. Then the equations for fluctuations 
AM2a;,y ~ exp(— zciji) reduce to 



iujM' 



2z 



AM. 



2y 



^2 , 
— AM. 



2y 



a 

iujM2z 
a 



AM2X = 0, (57) 



where two characteristic frequencies are introduced: 

ni = -i[HM2z+Ha+^nM2 
^2^1 (hM2z + i/a) - iKW. (58) 

We need now to linearize boundary conditions and use 
them together with equations (|57|) . Parameters k and p 



(|47|I " H49(I should be taken for the static state chosen in 
the linear approximation. It means we need to rewrite 
general expressions H49f) - (|51|l in a more specific form, 
which appears valid for any sign of current: 



k = 



] 



aM2 1 + h' e ' 



P 



llBlOi2T2Ql L 



(59) 



(60) 



Linearized boundary conditions after the transforma- 
tions mentioned acquire the following form: 



aAM2, 



dx 



dAM- 



2y 



dx 



x—Xf 



x—Xji 



kAIVhyiXp) - pM2zAM2x{\f), (61) 



-kAM2x{XF)-pM2zAM2y{XF), (62) 



dAM2x 



dx 



9AM. 



■2v 



x—L~e 



dx 



= 0. 



(63) 



x—L — e 



Let us substitute AM2x,y ~ eyi]i{iqx') into equa- 
tions H57|l and define x' = x ~ Xp- Then we find four 
possible characteristic wave numbers: 



4 



1 

2a 



(64) 



Then, the general solution of our linearized equa- 
tions lf57|) may be presented as follows 

AM2x{x') = Acosq-x' + Bs'mq-x' 
+C cos q+x' + D sin q+x\ 
AM2y (x') =^{A cos q_x' + B sin q_ x' ) 

-|-^^"^(Ccosq+x' + Dsing+x'), (65) 

where we took Re q± > and introduced a "distribution 
coefficient" 



i = 2iujM2z 



f^2 - Hi + v/(f^i-f^2)2+4a;2 



(66) 



Constants A, B, C, D should be found now from the 
boundary conditions (|61|l - (|63|l . Substituting H66() into 
the conditions and equating the determinant to zero, we 
find the dispersion relation for the fluctuations in the 
most general form: 

(^q^ tSLnq^L + pM2z^ (^q+tanq+L + pM2z 
2^fc 



+k 



[q- tan q^L — q+ tan q+L) = 0. (67) 



We may compare Eq. (I67|l with before derived disper- 
sion relation in Refs. Note, parameter 7^ 1. It 
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can be seen directly from the definitions H66() and (|58|l . 
Taking this into account, we see the relations mentioned 
may be exactly reduced to each other but for the case 
k = only, that is for the case of the LSI effective field 
mechanism is dominated. It appears due to rather crude 
approximation used in Refs. 718- when deriving the disper- 
sion relation. Contrary to this, our last expression Ht)7|) 
represents more correct and complete dispersion relation, 
which describes joint action of TST torque and LSU ef- 
fective field mechanisms. 



IX. INSTABILITY OF FLUCTUATIONS 

We need further to solve the dispersion relation ll^ 
and find eigenfrequencies uj for fluctuations. In particu- 
lar, it allows us to determine the current driven instabil- 
ity condition, that is the condition when Imw > 0. Gen- 
erally speaking, the transcendental equation Ht)7|) may be 
solved numerically for a number of experimentally given 
junction parameters. However, the actual task now is to 
develop a general representation about the switching pro- 
cess rather than to get specific numerical results. To fulfil 
the task, we may try to solve the equation analytically 
but only for a special case of nearly uniform fluctuations 
which satisfy the equality tang__L sa q^L. 

It is convenient to use Eqs. H64I) . (|66|) and rewrite 
Eq. Hfi7|) in the form 



pM2z ' 



-I- 



2a 

T 

2iujkM2 



a{ql 



1 - 



q+ t&nq+L + pA'hz 
qlL+pM2z 



q+ tan q+L + plVhz , 



(68) 



If current is absent (that is k 
Eq. ^ 



: p = 0), we obtain from 



\ 



4:TTjM2 
1 + 



hM- 



2z 



4^ 



(69) 



where we remain only positive frequency Rew > and 
introduced dimensionless fields: h = H/4:TtM2 and ha = 
Ha/4:TTM2. According to typical values of parameters, 
we have H ~ 10-100 Oe, Ha ~ 100 Oe, Ma ~ 10^ G 
and therefore h, ha < 0.01. We neglect small corrections 
originated due to higher powers of h, ha in Eq. H69|) . As 
it is seen from Eq. H69|l . the fluctuations without current 
are, of course, stable due to the Gilbert damping. 

Return now to the situation when the current is turned 
on (j ^ 0) but is not too large in magnitude to allow pre- 
serve in Eq. H68|) only linear in current summands. The 
justification of the supposition will be considered later. 



at the end of the section. If we accept now that the sup- 
position is valid, it becomes possible to insert q^ — 
into the square brackets of Eq. H68|l . It means, in partic- 
ular, we may lay (fii — fia)^ + = Qi + fl2 in the 
definition H64|) of and therefore use the expression 



^1+^2 



(70) 



to substitute it to the right hand side of Eq. (|68|l . More- 
over, we should neglect summands proportional to fc^ and 
kp too. Formally, these simplifications become valid if the 
following conditions are fulfilled: 



a\p\ 



L\ni + n2\ 



« 1, 



pq+ 



p_ 



< 1. 



(71) 



Then we resolve Eq. (|68|) with respect to square root 
and square the result. Some summands appear contain- 
ing second (or higher) power of the current {k^, kp and 
so on). These summands are small also and should be 
neglected due to the conditions H71|) . After we make the 
simplifications and substitute (|70() , we obtain the follow- 
ing square equation for the eigenfrequency uj: 



Re(rJi -t-r^a) - 



2apM2z 2akM2z 



Re (rJi^ia) - ^^^Re {fli + ^2) 



L 

= 0. 



kL 



(72) 



We should solve the equation (|72|) using the definitions 
of the parameters k and p given by Eqs. (|59|l and H6()|l . 
Then we obtain finally 



_^ _ M2z{j/e) 



.27k: 



K 
Kq 



ijo/\e\) 

M2zh _ M2z\j 
ha 

M2z{j/e) 




Uo/\e\) 

where we used the following simplifying assumption 

|fc| 



k\p\ 



» 1, 



(73) 



(74) 



which is really true for many actual cases of small enough 
Gilbert parameter k. Formula H73|) reduces to (|69|l in the 
limit j — > 0. The following additional system character- 
istics are introduced in Eq. (|73|l : 
characteristic damping 



2hav{i' + 1) 



7a2T2M2 Xvi{v-l)-2bvM'^ 



2z 



(75) 
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and 



characteristic current density 



{jo/\e\) = 



AnAhbhaXiiy + 1)^ 



MBa2T2Qi Ai^i(i/ - 1) - 2bvM2 



(76) 



Here the matching parameters vi and v are exactly 
the same as they were defined after formulae (|20|l . 1)21(1 
and (EEll- We may substitute the following typical val- 
ues: 1^2 ~ 1, f^i S> 1 and v ^ 1, \ ^ 0.1, Z ~ 17nm, 
aa ~ 2 X 10-^, ra - 3 x lO^i^s, ^ ^ 0.3. Then we 
get the estimations: |ko| ^ 3x10^^, |jo| ^ lxlO''A/cm2. 

Expression H73|l allows write the following two condi- 
tions of instability (Imo; > 0): 



- I?! 



Jo 



M2zh 

ha 



A^2.-^ > — . 
Jo/|e| Ko 



(77) 



(78) 



The fulfillment of any of the conditions (|77|l or H78|) 
or both of them together is sufficient for instability. The 
condition H77|l originates from p parameter and therefore 
refers to LSI mechanism. The condition l|78|l originates 
from k parameter and refers to TST mechanism. The 
threshold current jth has a minimal module, which is 
necessary to reach the instability according to Eqs. 177|) 
and ITHJl, that is 



jth = bo I 



M2,h 



ha 



■ (79) 



As it is seen, the threshold for LSI mechanism does 
not depend on the Gilbert constant n, while the TST 
threshold is directly proportional to it. Because of ex- 
perimentally revealed damping at room temperature"^^ 
K w (2 5) X 10^2 corresponds to the above given esti- 
mation K ^ \kq\ ~ 3 X IQ-^, we may expect any of the 
mechanism, in principle, is responsible for the instability. 
Note, the parameters jo and kq, according to the defini- 
tions H75II and H7t)|l . may change their signs and become 
negative. But it may occur only simultaneously for Jq 
and kq. Therefore, the TST condition (|78|l does not de- 
pend on the sign at all, but to provide the LSI instability 

condition be fulfilled we should change the sign of Al2z 
only. 

Let us take for further estimations: 



1 + 



M2zh 



ha 



< 1 or 



l«:o| 



K 

l«o| 



(80) 



and moreover j/ ^ 1, 1/2 ^ 1 and all the estimations from 
the text are valid. Then, we substitute these parameters 
into conditions H71|l and get 

a\p\ 



10-2, 




- lO-^ 


p_ 


- 10-2 




pq+ 




Q+ 





We see, therefore, the conditions H71() are good satis- 
fied for the threshold current and even above it. The last 
estimations justify our above-stated analytic approach to 
the solution of the dispersion relation Ht)7|) . It is interest- 
ing to stress, our eigenfrequency lu (|73() appears essen- 
tially nonlinear function of the current j /jo in spite of 
the linearization conditions H71|) . 



X. STATE DIAGRAM AND HYSTERESIS 



We discuss now in detail the instability conditions 1(77(1 
and ((78(1 . Let us draw the corresponding threshold curves 
in the current to damping plane — the "state diagram" . 
The result is presented qualitatively in Fig. |3| We con- 
fined ourselves only the case of positive parameters jo > 
and kq > and built up the boundaries of the regions in- 
dicated. The condition ((77(1 can be satisfied for M2Z = 1 
only. It means the LSI mechanism can excite the insta- 
bility of AP state only to switch it into P state. That 
is true whatever sign is taken for the current. Another 
situation appears in connection with the condition 1178(1 . 

For forward current j/e > the only value M2Z = 1 is 
acceptable, but for backward current j/e < it is neces- 
sary to have M2Z = — 1- It means the TST mechanism 
excites the instability of different states depending on the 
sign of the current: forward current excites AP state, but 
backward current excites P state. All these possibilities 
are seen from Fig. O 

Two notes should be made here. First, let us go along 
the line oia having AP state. In the point a the state 
becomes unstable and will be switched to the stable P 
state. We may return then to the point b, where the P 
state becomes unstable and will be switched to the initial 
AP state, which now is stable, and so on. Any switch- 
ing process leads to change of the junction resistance, 
the higher resistance being correspond to AP state and 
the lower resistance to P state. Such a behavior is illus- 
trated in Fig. ^ and similar observations were revealed 
experimentally in many papers. 

The second note concerns the intersections of our line 
oifl in Fig. 131 or some other lines (02C, for example) with 
the threshold curves in the points 6', d' and so on. In 
the points all the initial states, AP and P, become un- 
stable. Therefore, no switching can appear and instead 
some new time dependent oscillatory state, apparently, 
should form. This is a very interesting possibility, which 
correlates with the well known experimental observations 
of a microwave emission from the junctions above the in- 
stability threshold (see, e.g., Ref . [l31 . We plan to inves- 
tigate the possibility in more detail elsewhere. 
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FIG. 3: State diagram. Threshold current versus damping 
(qualitatively) at vi{v — 1)A ^ 2bv, vi 2> 1, 1^2 = 1. Heavy 
lines show boundaries of the regions. The following nota- 
tions are introduced: AP/LSI means that AP initial state 
is unstable due to LSI mechanism; P/TST means that P 
initial state is unstable due to TST mechanism, and so on. 
The region around the line j = without any designations 
corresponds to stability. Points of switching: a — switch- 
ing AP ^ P, 6 — switching P AP, c and d — switching 
AP P. Points of oscillations: b' and d! . 



R 



jle 



FIG. 4: Hysteretic behavior of junction resistance R: Rap, 
Rp are resistances of AP and P states, respectively. 



XL SPECTRUM SOFTENING AND 
INCREMENT 

We discuss now he question: how spectrum (Rew) and 
increment (Im w) behave when the current comes near the 
threshold of instability. To answer, we consider Eq. H73|) 
for eigenfrequency and calculate numerically the quan- 
tities of interest for some typical example. Figs. |3HZ| 
demonstrate the most significant results. 

According to Fig. |S1 a normalized spectrum Rcw/wq 




FIG. 5: (Color online) Current dependence of eigenfrequency 
Retj/ojo (curve 1) and decrement Imtj/cjo (curve 2). Damp- 
ing constant /v/kq = 0.5. 



falls to zero, that is softens, when the current rises near- 
ing to the LSI threshold. As it was remarked in Ref. 0, 
the softening appears due to LSI driven reorientation 
phase transition at the threshold. Indeed, such a soft- 
ening is the general property of similar phase transitions. 
Contrary to this, the TST has no relation to any phase 
transition^ and is simply a mechanism to deliver energy 
into the magnetic lattice, correspondingly, increment 
Imw/wo changes its sign and becomes positive not at 
the LSI threshold but a little below if damping is small 
enough, namely, for k/kq = 0-5 in our example. It is 
significant to stress the increment rises faster when cur- 
rent increases in the region where the LSI mechanism 
dominates. If we estimate ojq ^ 10^" s~^, then the in- 
crement may rise from approximately Imw ^ 10^ s""'^ for 
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10 , 



for LSI 



TST mechanism dominated to Imw 
mechanism dominated. 

Fig. |B| shows the spectrum for different values of an 
external magnetic field. Eigenfrequencies rise with in- 
creasing the field. The LSI threshold shifts toward large 
currents. But the fact of spectrum softening, of course, 
remains. 

Fig. El shows the increment for different values of a 
damping constant and in absence of external magnetic 
field. As the damping increases, the increment at a given 
current decreases for any mechanisms (TST and LSI) but 
LSI threshold is fixed, that is it does not depend on the 
damping. The TST threshold rises with damping and 
reaches the LSI threshold at the end. 



XII. CONCLUSION 

A phenomenological approach is developed in the the- 
ory of spin-valve type ferromagnetic junctions. New vec- 
tor boundary conditions are derived, which represent the 
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FIG. 6: (Color online) Softening of eigenfrequency for differ- 
ent external fields: 1 — h/ha = 0, 2 — h/ha ~ 0.5, 3 — 
h/ha = 1.0. Parameters: ^ 1, b ^ 1, K,o/{2\/hI) = 0.1. 



Im oj/to^ 



1.2- 1 




FIG. 7: (Color online) Increments for different damping: 1 — 
k/ko = 0.5, 2 — k/kq = 1.0, 3 — k/ko = 2.0. Parameters: 
> 1, b ~ 1, ft = 0, Ko/(2V7t") = 0.1. 



continuity conditions for different spin fluxes: longitudi- 
nal spin flux of mobile electrons and transverse spin fluxes 
of the magnetic lattice including the so called "effective 
field transverse spin flux" . Spin transfer torque is con- 
sidered as a boundary condition. It describes the result 
of abrupt transformation of the mobile electron trans- 
verse spin fiux to the spin fluxes of the lattice. Joint 
action of two electric current effects is investigated: the 
nonequilibrium longitudinal spin injection effective field 
and the transverse spin transfer surface torque. General 
macroscopic dynamic equations are linearized for fluctu- 
ations and solved using the boundary conditions men- 
tioned. Dispersion relation is derived and solved. 



The calculation model is extended significantly to 
make it adequate the experimental situation as far as it 
is possible. In particular, the mobile electrons are consid- 
ered as having non-pinned spins in all the contacting lay- 
ers. Lattice magnetization is assumed completely pinned 
in one of the magnetic layers and free in the other mag- 
netic layer. Theory is developed for any current direction 
in the spin- valve junction (forward or backward). 

Some critical value kq of the well known Gilbert damp- 
ing constant k is introduced for the first time. Spin trans- 
fer torque gives the dominant contribution to instability 
threshold for small enough constants k < kq, while the 
spin injection contribution dominates for k > kq. Typi- 
cal estimation gives kq ^ 3 x 10~^ and it approximately 
corresponds to the experimentally found Gilbert constant 
values for ferromagnetic metals at room temperature. It 
means both the two mechanisms may be responsible for 
the measured thresholds. 

The mechanism of magnetic switching by the back- 
ward current and hysteretic behavior is proposed for the 
first time. As it appears, a fine interplay of spin transfer 
torque and spin injection is necessary to provide the hys- 
teresis. The peculiarity of the process may be understood 
as follows. Spin injection effective field is localized in the 
free layer near the boundary between pinned and free 
layers. Spins of mobile electrons in the backward current 
have a component transverse to the effective field and 
therefore perform a precession in the field. This preces- 
sion leads to torque generation and switching in complete 
accordance with the transverse spin transfer mechanism. 
The consequent realization of forward and such type of 
backward switching leads to hysteretic behavior shown in 
Fig. 21 The existence of the hysteresis is consistent with 
experimental data. 

The state diagram was built up that divided the 
current-damping plane (Fig. O into regions of fluctu- 
ation stability and regions of various type instabilities. 
For example, for forward currents the AP initial states 
may become unstable and P states remain stable. In the 
regions the switching processes (i.e., AP P) become 
possible. This diagram revealed for the first time the 
existence of another type regions for backward currents 
where AP and P initial states become unstable simul- 
taneously. This last type of regions may correspond to 
formation of time dependent (oscillatory) states as a re- 
sult of the instability development. 

The effect of softening was predicted for the layer spin 
wave eigenfrequencies when current rises and comes to 
the instability threshold (Figs. [S] and 0). This effect is 
interesting because it distinguishes directly between TST 
and LSI mechanisms, since only the LSI mechanism leads 
to the softening. 

The estimations were made of instability increments as 
functions of the current. It was shown for the first time 
that LSI mechanism leads to much larger increments in 
comparison with TST mechanism. For example, if we 
take spin wave eigenfrequency Rew 10^° s~^, then we 
get increment Imw ~ 10^ s^"'^ for TST mechanism dom- 
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inated and Imo; 



10 



10 , 



for LSI mechanism domi- 



nated. The increment rises for any mechanism if current 
increases and damping decreases. 
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APPENDIX: TRANSVERSE SPIN TRANSFER 
TORQUE CONSIDERED AS A LOCALIZED OR 
AS A UNIFORMLY SPREAD 

We try now to show the connections between two mod- 
els discussed till now: traditional model of uniformly 
sprea d TST torque and recently introduced (Refs. 011, 
IgllOl and this paper) localized TST torque. 

We start from the linearized equations and bound- 
ary conditions for fluctuations (|57|l and (|61|l - (|63|l . The 
boundary conditions (|61|l and H62|l may be inserted into 
Eq. H57() using a singular J-function, namely, 
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-kAAhx - pM2zAM2v\ S{x - Xp). (A.l) 



To check the correspondence of the Eqs. ljA.l|) 
and ((23, (EB, (E2Il, let us integrate (|A.1|1 over the in- 
terval [0, Ai? + 0]. Then, we have for the derivative 
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(A.2) 
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and analogous expression for the other derivative, which 
may be got by a replacement x y, where we accept 
free spins at the plane x = 0, that is 



dAM2x 



dx 



dAM. 
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x=0 



dx 



= 0. 



(A.3) 



x=0 



Moreover, we should neglect all the terms, which are 
proportional to small length — > 0. Then we get imme- 
diately the boundary conditions in the form H61|l . (|62|) . 



Therefore, two completely equivalent approaches may be 
used: either we solve our previous equations l|57|) satisfy- 
ing nonuniform system of boundary conditions (|61|) - (|63|l . 
or we solve the singular equations l|A.l|) with uniform 
boundary conditions and HA.3|I . 

The second approach appears convenient to pass to 
uniform model. To do this, we should spread uniformly 
through the whole layer 2 thickness L the localized exci- 
tation in ljA.l|) . that is we should take 



S{x - Xf) —. 



(A.4) 



After this is done, the equations IIA.l|) become spatially 
uniform and exactly coincide with the linearized uniform 
model. The last statement may be directly confirmed if 
we start from LLG equation with the well known addi- 
tional double vector product term introduced in the orig- 
inal worksiiS. The boundary conditions (|63|l and ljA.3|) 
are uniform also. It shows we can find a uniform solution 
for fluctuations. To simplify extremely the calculations, 
we suppose the thickness of the layer is sufficiently small 
and spin injection may be omitted. Then the following 
dispersion relation for eigenfrequency becomes valid in 
linear in current approximation: 



(1 + k^)llP' + inuj 



Re{n, + n2)-^M2z 



Re(17ir22) = 0. 



(A.5) 



This dispersion relation coincides with the correspond- 
ing relation of our localized torque model H72|l taken at 
p = 0. The condition of instability that follows from 
(|A.5p may be written in the form 



M2 



jth.± ^ 27r7M|Ki(l + I') 



(A.6) 



As it is seen, forward current {jth,±/s > 0) leads to in- 
stability of AP states only (because of to satisfy (|A.6p we 

should take M2Z — 1) and, vise versa, backward current 
{jth,±/e < 0) leads to instability of P states (we should 

take M2Z = — !)• This conclusion is in agreement with 
our localized torque theory (see Section^of this paper). 
Numerical estimation gives \jth,±/e\ ~ 10^ A/cm^. 

However, we should stress here the correct understand- 
ing of the problem cannot be possible without including 
both the mechanisms (TST and LSI) acting simultane- 
ously. First, the LSI was included implicitly in the con- 
dition (jA.6p for backward current switching. Really, the 
existence of sd exchange effective field AHsd ^ d{x — XF) 
(see Eq. (|25ll ') is necessary to understand the backward 
transverse spin flux appearing in the layer 2 and transfor- 
mation of this flux into the lattice flux. The relation 144|) 
describes such a transformation and allows calculate the 
TST parameter k (|5^ . 

Moreover, the results of the paper show the LSI mech- 
anism manifests itself significantly and in different ways. 
For example, the instability threshold does not rise with 
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increasing the Gilbert damping constant k but reaches 
some Umiting value due to LSI mechanism (see Fig. O. 
Increment of instability may be more than order of mag- 
nitude larger due to LSI (see Fig. (SJ. Interplay of TST 
and LSI mechanisms makes it possible to provide the in- 



stability simultaneously of P and AP states (Fig.lSJ. It 
means, no switching but rather new time dependent final 
nonlinear state may appear as a result of the instability 
development. 
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tween mobile electron spins and bound electrons providing 
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ferent ways of the action. It will be treated further in more 
detail. 
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media approximation. 
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might arise due to spatial distributions only: random or 
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